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We present a fit to measured moments of inclusive distributions in B — > X c lv and B — > X s 7 
decays to extract values for the CKM matrix element |V c t|, the b- and c- quark masses, and higher 
order parameters that appear in the Heavy Quark Expansion. The fit is carried out using theoretical 
calculations in the kinetic scheme and includes moment measurements of the BABAR, Belle, CDF, 
CLEO and DELPHI collaborations for which correlation matrices have been published. We find 
| Vci,| = (41.96 ± 0.23 cxp ± 0.35 H qe ± 0.59r SL ) x !0~ 3 and m b = 4.590 ± 0.025 cxp ± 0.030 H qe GeV 
where the errors are experimental and theoretical respectively. We also derive values for the heavy 
quark distribution function parameters mj and in different theoretical schemes that can be used 
as input for the determination of \V u t\- 

PACS numbers: 12.15.Ff, 12.15.Hh, 12.39.Hg, 13.30.Ce 
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I. INTRODUCTION 

In the past few years tremendous progress has been 
made in the description of semileptonic and radiative 
B decays using the framework of Heavy Quark Expan- 
sions (HQEs). Calculations for the semileptonic decay 
width as well as for moments of inclusive observables 
with restrictions on the phase space are now available 
in different schemes through order l/m\ and oqPq 
Q 13- S D- IS ■ At the same time many new experimen- 
tal measurements of moments of the hadronic mass and 
lepton energy distribution in B — > X c lv as well as the 
photon energy spectrum in B — > X s j decays have been 
carried out by several experiments. Generally, the results 
agree very well between experiments. In addition, the 
theoretical calculations describe the measured data well 
establishing this framework for treating semileptonic 
and radiative B decays 0, IS El 

In this document we will present the results of a com- 
bined fit to measured moments for which correlation 
matrices are published. These include moments of the 
hadronic mass distribution (M^) and moments of the 
lepton energy spectrum (Ef) in inclusive B —>■ X C £D de- 
cays as well as moments of the photon energy spectrum 
(E") in inclusive B — > X s j decays for different mini- 
mum lepton and photon energies E cut - The HQEs for 
the moments depend on the b- and c- quark masses and 
several non-perturbativc parameters which therefore can 
be determined from a fit of the theoretical expressions 
to the experimental moment measurements. Among the 
measurements we have excluded those for which there 
are ongoing discussions within the theoretical commu- 
nity regarding the associated theoretical uncertainties. 
These are in particular the non- integer moments (Mj) 
and (Mjf ) of the hadronic mass distribution in B — > X C £D 
decays. In addition, moments of the photon energy spec- 



trum in B — > X s j decays above E cut = 2.0 GeV have 
been excluded since there the standard local Operator 
Product Expansion (OPE) is no longer believed to be 
under theoretical control. Furthermore, individual mo- 
ment measurements have been discarded if the covariancc 
matrix cannot be inverted due to the large correlations 
between the measurements. 



II. HEAVY QUARK EXPANSIONS IN THE 
KINETIC SCHEME 

In this analysis we make use of HQEs that express 
the semileptonic decay width TgL as well as moments 
of the lepton energy and hadron mass distribution in 
B — > XJLv decays and those of the photon energy spec- 
trum in B — > X s "f decays in terms of the running kinetic 
quark masses m&(/x) and m c (/i). Non-pcrturbative effects 
are introduced in this formalism via heavy quark oper- 
ators. The leading power corrections arise at 0(1/ ml) 
and are controlled by the two expectation values 
and /.Iq (fi) of the kinetic and chromomagnctic dimension- 
five operators. At 0(l/m|) two additional expectation 
values p%(n) and p^sil 1 ) °f tne Darwin and spin-orbital 
(LS) dimension-six operators complete the set of non- 
perturbative corrections. Together with the two running 
quark masses the HQE in the kinetic scheme includes six 
free parameter through 0(l/m^): 

• Leading Order Parameters 

— m(,(/i) — > b-quark mass 

— m c ([i) — » c-quark mass 

• Leading Non-Perturbative Corrections - 0(1/ ml) 

— fJ^ifM) — > 'kinetic expectation value' 

— Mg(/i) — > 'chromomagnctic expectation value' 
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• Higher Order Non-Perturbative Corrections 
0(l/mg) 



Pniv) - 
Pls(p) 



'Darwin term' 
'spin-orbital term' 



At any given value of the Wilsonian factorization scale 
fj, all of the above mentioned HQE parameters represent 
well defined physical quantities which have to be deter- 
mined from measurements. The scale fi separates 'short- 
distance' effects from 'long-distance' effects and, there- 
fore, rriQ(fi) can be understood as a short-distance mass 
of perturbation theory that excludes soft gluon interac- 
tions 0. It is important to note that Ref. [j| contains 
a translation to full order a 2 and third-order BLM cor- 
rections of the running short-distance mass 771q(/x) into 
the well known running MS mass ^(wq). Hence the 
kinetic masses mb(p) and m c (/i) can be compared with 
other established mass definitions in QCD such as the 
MS mass. This fact is important for the comparison with 
other QCD calculations beyond semileptonic or rare B- 
decays. In order to minimize the influence of radiative 
corrections and to insure that the kinetic c-quark mass 
m c (p) has a well defined physical meaning the separation 
scale is set to be // = 1 GeV. A detailed discussion of the 
justification for this choice can be found in Ref. [lj}. ^ n 
the following, all heavy quark parameter values are pre- 
sented for fx — 1 GeV (e.g. m;,(l GeV) — ♦ mb). The ana- 
lytical expression for the semileptonic width of B — > X C £D 
decays through 0(1/ ml) Q is given by 

T8h(B->Xctv) = ^4\V cb \ 2 (l + A ew )A pert (r,p) 




(1) 



d(r)^ + 0(l/mt) 
m 3 



where r = m 2 /m 2 explicitly contains the c-quark mass 
m c . The tree level phase space factor z$ is defined 
through: 

z (r) = 1 - 8r + 8r 3 - r 4 - 12r 2 ln(r) , 

and the expression d(r) is given by: 

, / s 34 32 , 32 , 10- 
d(r) = 81n(r) + — - — r - 8r 2 + —r 3 - —r 4 . 

The electroweak corrections for Eq.n arc estimated to be 
approximately 



l+A e 



a M z , 
I + -In— ) w 1.014 

7r mb 



and the pcrturbative contributions, which have been 
calculated to all orders in BLM corrections and to second 
order in non-BLM corrections, are for a reasonable set 



of HQE parameters estimated to be A pert ~ 0.908. It 
should be noted that Eq. ^ is n °t a HQE in powers 
of ljrtib- Since the most relevant scale for the b — > c 
transition is the energy release mj, — m c of the decay 
rather than the b-quark mass mb, this expansion is 
carried out in powers of \/(mb — m c ). On the other 
hand, due to the low mass of the s— and u— quark, HQE 
calculations for the b — > s or b — > u transition can be 
considered as expansions in powers of 1/mfc. The full a 2 . 
corrections to the photon spectrum have been computed 
recently ^ s follows from this result, the effects 

of omitting the non-BLM corrections are small and 
fully covered by the assumed theoretical uncertainties 
we quote for the first and second photon energy moments. 

For the practical use Eq. ^can be transformed into: 



WrJ 



BRciv 



TB I 



0.105 - 0.0018, 
x [1 + 0.30 (a s (m&) -0.22)] (2) 
x [1-0.66 (m;, -4.6 GeV) + 0.39 (m c - 1.15 GeV) 
+ 0.013 (/4 -0.4 GeV 2 ) +0.09 (po 3 -0.1 GeV 3 ) 
+ 0.05 (fi% -0.35 GeV 2 ) - 0.01 (pl s +0.15 GeV 3 )] , 

where \V® b \ — 0.0417 and tb represents the average life- 
time of neutral and charged B mesons. We use tb = 
1.585 ± 0.007 ps based on Ref. [13, assuming equal pro- 
duction of charged and neutral B mesons. The theoret- 
ical uncertainty due to the the limited accuracy of this 
HQE is denoted by <5 th . In Ref. its value is quoted 
with (5 t h = 0.015, with contributions from four different 
sources. Since then a more elaborate study on the influ- 
ence of 'intrinsic charm' has been carried out redu cing 
the associated uncertainty by roughly a factor two 
leaving an overall uncertainty of 5 t h = 0.014. More de- 
tails about this theoretical uncertainty can be found in 
Section ITTTCl 

It should be noted that we do not extract the Dar- 
win term p|)(l GeV) directly from our fit to the HQEs 
but rather the 'pole-type' Darwin expectation value pb 3 ■ 
The two parameters are closely related in our framework 
via 



PD 3 



pfj(lGeV) -0.1 GeV 3 



Relations similar to Eq. have been calculated for in- 
clusive obscrvablcs in semileptonic and radiative B de- 
cays. These are in particular hadron mass (M^) and 
lepton energy moments (Eg) to order n as well as mo- 
ments of the photon energy spectrum {E™): 



(Ml) 
(Ef) 



( M x) (mb, m c , /4, /4; Pd,Pls, a s 
(Ee) (m b , m c , /i 2 , /1q, p 3 D , p 3 LS , a s ) 



Since every moment calculation has a different depen- 
dence on the heavy quark parameters a simultaneous fit 
allows for the extraction of all these parameters. For this 
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TABLE I: Summary of moment measurements used in the combined fit. n indicates the order of the (central) moment 
measurement of observable (M^)_E cut , (iJ™)_B cut and {E")E CU f E cut indicates measurements with the corresponding minimum 
lepton momenta and photon energies in GeV. 



Experiment 


Hadron Moments (Mx}b oM 


Lepton Moments (_E").E cut 


Photon Moments (E") Ecllt 


d a d ad 

[iftiJJ 


71= z _C/ cut =U.y,l.U,l.l,1.2,l.o,1.4,1.0 
n=4 £ cut =0.9,1.0,l. 1,1.2,1.3,1.4,1.5 


n=U -£/ cu t=U.O,1.2,1.0 

n=l £ cu t=0.6,0.8,1.0,1.2,1.5 
n=2 £ cut =0.6,1.0,1.5 

„ Q T? f\ Q 1 

n— ,i iicut— U.o,l.z 


n=l _£/ cut =1.9,2.U 
n=2 £ cut =1.9 a 


Belle 
[18, Wj 






n=l £ cut =1.8,1.9 
n=2 £ cut =1.8,2.0 


CDF 
[20] 


n=2 B cu t=0.7 
n=4 E cut =0.7 






CLEO 
[21, 22} 


n=2 £ cut =1.0,1.5 
n=4 £ cut =1.0,1.5 




n=l £ cu t=2.0 


DELPHI 

w 


n=2 E cut =0.0 
n=4 E cut =0.0 
n=6 £ C ut=0.0 


n=l £ cut =0.0 
n=2 £ cu t=0.0 
n=3 £ cu t=0.0 




HFAG [23J 




n=0 £ cut =0.6 





a A total of six photon moments from Refs. I Hi and are used. 



it is important to use as many moment measurements as 
possible in order to overconstrain the extraction of the 
heavy quark parameters and to establish the validity of 
the expansions. A much more detailed description of the 
theoretical framework used for this analysis can be found 
in Refs. HHH- 



III. FIT TO MOMENT MEASUREMENTS 

In the following sections we list the currently available 
experimental moment measurements and indicate which 
are used in the combined fit presented here. Furthermore, 
we outline the fit procedure and summarize the results 
obtained from different fit scenarios and various cross 
checks. 



A. Experimental Input 

All results are based on the following set of moment 
measurements which are also summarised in Tabled Ad- 
ditional measurements for which correlation matrices are 
not available and thus cannot be used in the presented 
fit are listed in parentheses. 

• BABAR 

Hadron mass 0] and lepton energy moments 0] 
from B — > X C £D decays measured as a function of 
the minimum lepton energy E cut - The lepton mo- 
ments used here differ slightly from those in the 
BABAR publication 0. They have been updated 
by taking into account the recent improved mea- 
surements of the D s and B branching fractions 
{upper-vertex charm) that impact the background 



subtraction. Moments of the photon energy spec- 
trum in B — > Xsj decays as a function of the min- 
imum photon energy E cut from two independent 
analyses [11,113 

• Belle 

First and second moment of the photon energy 
spectrum as a function of the minimum photon en- 
ergy E cut MM 

(Measurements of hadron mass and lepton energy 
moments as functions of the lower lepton energy 
exist 0, H3 but are excluded from the current fit 
as correlation matrices are only available for the 
statistical errors.) 

• CDF 

Hadron moment measurements with a minimum 
lepton energy of B cut = 0.7 GeV HJ. 



• CLEO 

Hadron moment measurements as a function of the 
minimum lepton energy |2lj| . 

First (and second) moment of the photon energy 
spectrum at B cut = 2.0 GeV M 
(The measurement of lepton energy moments as a 
function of E cut j2|| is not given with the full co- 
variance matrix and thus has not been included in 
the fit H3-) 



• DELPHI 

Lepton energy and hadron mass moment measure- 
ments with no restriction on the lepton energy 
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Fit Procedure 



A x 2 minimization technique is used to determine the 
HQE predictions in the kinetic scheme from a fit to the 
data: 

X 2 = (Mcxp - Mrqe) T C^I (M cxp - A/hqe) (3) 



exp represents all moment measurements 



where M, 

included in the fit and Mhqe stands for their HQE 
prediction defined by the heavy quark parameters. 
Ctot = C cxp + Cthco is the sum of the experimental and 
theoretical covariance matrices. The construction of the 
theoretical covariance matrix Cthco is discussed in detail 
in Section 

In order to extract the semileptonic branching fraction 
for B — > X c lv events, BR c gp, from the fit, the mea- 
surements of the truncated branching fractions (E^) at 
different cutoffs in the lepton energy E cut are also used 
as experimental input. HQE predictions of the relative 
decay fraction for a given cutoff -E cu t; 



R(Ecut) 



Jo 



dE L 



dEr 



(4) 



can be used to extrapolate the measurement of the trun- 
cated branching fractions to zero cutoff: 



BR C 



{ED 



R(E cu t) 



(5) 



In addition to the BABAR measurement we also 
include the HFAG average at E cut = 0.6 GcV of 
(10.29 ± 0.18)% In order to utilize more than only 
one of these truncated branching fractions, one has to 
include the total branching fraction BR C £ P as a free 
parameter in the fit. Together with the input of the 
averaged B meson lifetime tb = 1.585 ± 0.007 ps this 
can be used to calculate the semileptonic width Tg^ 
as part of the \ 2 minimization. The HQE for is 
directly related to \V c b\ 2 (see Eq. 0) and introducing 
\V c b\ as a free parameter therefore has the advantage 
of determining the error on this quantity directly 
from the global fit. With this approach also potential 
non-Gaussian errors (e.g. asymmetric errors) of the fit 
parameters are properly propagated into the error on 

\Vcb\- 

The fit to the moment measurements is carried out 
using the HQE calculations in the kinetic scheme pre- 
sented in Refs. @>B, including E cut dependent pertur- 
bative corrections to the hadron moments 0, l28l l29l |. 
However, rather than using linearized tables to deter- 
mine Mhqe (as was done previously in Ref. Q) , we 
obtain the prediction for every single moment from an 
analytical calculation [3(|. This not only allows us to 
study the scale dependence of the kinetic scheme in de- 
tail but also provides more accurate predictions for the 



individual moments. Since some criticism has recently 
been raised concerning the quality of the theoretical ex- 
pansion for the non-integer hadron moments (Mx) and 
(M x ), we exclude these moments from the fit until the 
issue of the theoretical uncertainty of these moment pre- 
dictions is resolved. However, for comparison, we will 
always show their prediction based on the fit results and 
compare them with the corresponding measurements. 

As ixq and p^ s arc estimated from the B* — B mass 
splitting and heavy-quark sum rules, respectively, we im- 
pose Gaussian error constraints of p? G = 0.35 ± 0.07 GeV 2 
and p\ s = — 0.15± 0.10 GcV 3 on these parameters as ad- 



vocated in Ref. 0] ■ 



C. Theoretical Error Estimates 

Since this analysis targets a measurement of the CKM 
matrix element \V c b\ with a relative error below the 2% 
level it is of vital importance to take theoretical uncer- 
tainties into account, as currently estimated. The HQEs 
for the moments have theoretical uncertainties due to cer- 
tain limitations in the accuracy of the calculations and 
certain approximations. As Eq. [21 illustrates, the heavy 
quark parameters extracted from the simultaneous fit to 
the moments are used for a residual correction to \V® b \. 
It is therefore important to include their theoretical un- 
certainties in the total covariance matrix of the \ 2 fit to 
achieve realistic error estimates for them and for | V c b \ ■ 

Since the HQE for Tsl is the best known expansion 
in this framework also its error estimates are the most 
advanced. In Eq. [21 the HQE for \V c b\ is explicitly 
quoted with a theoretical uncertainty S t h to reflect 
the limited accuracy of this theoretical expression. 
Adding the individual uncertainties quoted in Ref. 
in quadrature and accounting for the more advanced 
estimates of the potential importance of 'intrinsic 
charm' of Ref. yields 8th = 1-4%. We quote this er- 
ror separately for \V c b\ and label it with a subscript Tsl'- 

Following the recipe quoted in Ref. 0] we estimate 
the uncertainties in the individual HQEs for the mo- 
ments M by conservatively varying the values of /i 2 
and /j,q by ±20% and those of the 0(l/m 3 ) operators 
pfr S and pfrj by ±30%. These variations are carried out 
around the expected theoretical values of p 2 . = 0.4 GeV 2 , 
Mg = 0.35 GeV 2 , p 3 D = 0.2 GeV 3 and p\ s = -0.15 GeV 3 . 
Perturbative uncertainties are addressed by varying a s 
by ±0.1 for hadron moments and by ±0.04 for lepton and 
photon moments around a central value of a s = 0.22. Fi- 
nally, we also vary the quark masses nib and m c by 20 
McV around central values of 4.6 GeV and 1.18 GcV. 

These uncertainties are propagated into an error on 
the individual moments AM using Gaussian error prop- 
agation 



(am) 2 = E(|f ) 2 4 



(6) 
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where £ = (mj, m c , /i^., /Xg, /r^,, pf^, a s ). All these vari- 
ations are considered uncorrelated for a given moment. 
The theoretical covariance matrix is then constructed by 
treating these errors as fully correlated for a given mo- 
ment with different E cut while they are treated as uncor- 
related between moments of different order. 

For the moments of the photon energy spectrum we 
include additional theoretical errors related to the ap- 
plied bias corrections as described in Ref. We follow 
the suggestion of the authors and take 30% of the abso- 
lute value of the particular bias as its uncertainty |3ll |. 
In addition we linearly add half the difference in the mo- 
ments derived from the two distribution function ansiitze 
as given in Ref. Q. These additional theoretical errors 
related to the photon energy moments are considered to 
be uncorrelated for moments with different E cut and of 
different order. 

Generally, the chosen approach for the evaluation of 
the theoretical uncertainties is conservative. As a re- 
sult of this the x 2 /Adof is very good as will be shown in 
the next section. It is interesting to note though, that 
a fit that neglects the theoretical errors leads to similar 
results. However, we consider the choice of theoretical 
errors as appropriate since we are trying to consistently 
extract \V c b\ and the heavy quark distribution function 
parameters from a fit to all moment measurements, and 
aim to arrive at a reliable error on \ V u b\ when using these 
results as input. The conservative approach is also re- 
flected in the fact that we exclude non-integer hadron 
moments and photon moments above i? C ut = 2.0 GcV 
from the fit. 



D. Fit Results 

In the following we present the results of a combined 
fit of the HQEs to all moment measurements listed in 
Table [I] along with their corresponding theoretical error 
estimates as defined in Section lill CI 

In order to asses the consistency of the moment 
measurements from the two different decay processes, 
B — > X c iv and B — ► X s ~/ , we also carry out sep- 
arate fits to B — > X C £D moments and to photon 
moments only. However, as the latter are not sensi- 
tive to all the heavy quark parameters, all but m,b and 
^ are fixed to the result obtained from the combined fit. 

A detailed comparison of the HQE predictions ob- 
tained from the combined fit and the moment measure- 
ments is shown in Figures 1 1131 for the hadron, lepton and 
photon moments, respectively. The yellow bands repre- 
sent the total experimental and theoretical fit uncertainty 
as obtained by converting the fit errors for each individ- 
ual HQE parameter into an error for the individual mo- 
ment. The green band indicates the experimental error 
only. These figures also show the measurements that are 
not included in the fit for the reasons described in Sec- 
tion IIIIA1 In particular the non-integer hadron moments 



(Mx) and (M x ) can therefore be directly compared with 
the corresponding fit prediction. It can be seen that all 
moment measurements agree with each other and that 
the fit is able to describe all the moment measurements 
of different order and from different B decay distribu- 
tions. 

Figure [3] shows a comparison of the fit prediction 
for the first and second moments of the photon energy 
spectrum for the standard OPE ansatz with the E^t 
dependent bias corrected OPE calculations of Ref. [j]. 
While it is expected that the pure local OPE approach 
will break down at higher values of E cut this is the first 
time that the accuracy of the experimental moment 
measurements is sufficient to demonstrate this effect. 
For i? cu t above 2.0 GeV it is clearly visible that the pure 
OPE ansatz fails to describe the data. 

The fit results are summarised in Table ITU where the 
separation of the errors into experimental and theoretical 
contributions was obtained from toy Monte Carlo exper- 
iments. The results are in good agreement with earlier 
determinations @ ; B but have improved accuracy. A 
comparison of results from the combined fit with those 
obtained from fits to B — > XJLv and B — > X s j moments 
only can be found in Figure 0] where the A% 2 = 1 con- 
tours for the fit results are shown in the (nib, \V c b\) and 
(mj„/jj) planes. It can be seen that the inclusion of the 
photon energy moments adds additional sensitivity to the 
b-quark mass m\y. 

Including the Belle measurements of the hadron mass 
and lepton energy moments with only their statistical 
correlations leads to very similar results with only small 
improvements in the errors for the heavy quark parame- 
ters. This is a consequence of the fit errors being dom- 
inated by the theoretical uncertainties as can be seen 
from Table [HJ For a future average based on the full 
covariance matrix it will be necessary to consider the 
strong correlation of the systematic errors of the lepton 
moments between experiments. In particular, a consis- 
tent background modeling and subtraction technique will 
be required. 

To ensure the stability of the fit procedure several 
cross checks have been carried out. For instance, the 
combined fit has been repeated without applying the 
theoretical constraints on /j,q and p LS . We also repeated 
the fit excluding hadron moments with units GeV 6 and 
lepton moments with units GeV 3 as these moments are 
believed to have large theoretical uncertainties (as can 
bee seen from Figures ^ and [2J . In addition, photon 
moments with i? cu t > 1.8 GeV were excluded as here 
the bias corrections become noticeable. Finally a fit 
neglecting all theoretical errors was performed, i.e. only 
the experimental covariance matrix was used. All these 
results agree well with each other and any variations 
are fully covered by the theoretical error estimates. In 
addition the scale dependence of the expressions for 
the moments was studied but was found to be small 
compared to the assigned theoretical uncertainties. 
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FIG. 1: Comparison of fit predictions and the hadron moment measurements: (a) (Mx), (b) (M x ), (c) (M x ), (d) {{M x — 
(M|)) 2 ) and (e) ((M x — {M x ))' i ). The yellow bands represent the total experimental and theoretical fit uncertainty as obtained 
by converting the fit errors of each individual HQE parameter into an error for the individual moment. The green band indicates 
the experimental uncertainty only. Solid markers are included in the fit while open markers are only overlaid for comparison. 
Moment measurements at different E cu t are highly correlated. 



In addition to the above we extract the difference in 
the quark masses as 

m h -m c = 3.446 ± 0.025 GeV . 

Comparing the extracted values of the quark masses 
nib and m c with other determinations is often convenient 
in the commonly used MS scheme. The translation be- 
tween the kinetic and MS masses to two loop accuracy 
and includingthe BLM part of the corrections was 
given in Ref. yj- This leads to 

m{m) = 4.20 ± 0.04 GeV 
m^(m^) = 1.24 ± 0.07 GeV 

These results agree well with the determination in the 
IS scheme P a , l32l and recent unquenched lattice calcula- 
tions |33l I34L l35j | . However, it has been accepted among 



theorists that the normalization scale of around 1.2 GeV 
in the MS scheme may be too low for a precision evalua- 
tion of masses, and higher-order perturbative corrections 
in m c (m c ) are too significant. As a result, an additional 
uncertainty in rn^(m^) of at least 50 MeV may have to 
be added associated with the definition of rn^(rrT^) itself. 
A larger normalization scale for the MS masses is gener- 
ally used. To address this we give here the value of m c 
normalized at a safer momentum scale 2.5 GeV as was 
advocated recently: 

rf^(2.5 GeV) = 1.072 ± 0.06 GeV . 

The theoretical uncertainty in this translation is small. 
It may also be convenient to have the ratio of the charm 
and the beauty quark masses in the MS scheme which is 
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• BABAR ■ BELLE ★ DELPHI * HFAG 




E cut (GeV) E cut (GeV) 3-2006 

8730A8 

FIG. 2: Comparison of fit predictions and measurements for the lepton moments: (a) BR, (b) {El}, (c) ((El — (El)) 2 } and 
(d) ((El — (El)) 3 )- The yellow bands represent the total experimental and theoretical fit uncertainty while the green band 
indicates the experimental uncertainty only. Solid markers are included in the fit while open markers are only overlaid for 
comparison. Moment measurements at different E cu t are highly correlated. 



normalization-scale independent : 

= 0.235 ±0.012. 

The uncertainty in this ratio is dominated by the fit error 
on m r . 
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• BABAR ▼ BABAR 

■ BELLE ▲ CLEO 



OPE OPE + Bias Correction 




1.6 1.8 2.0 2.2 1.6 1.8 2.0 2.2 

E.(GeV) E cut (GeV) 10-2005 

cut CUt» ' 8730A3 



FIG. 3: Comparison of fit predictions and measurements for the photon moments: (a) (E^) and (b) ((E~, — (E^)) 2 ). The bands 
in the figures on the left show the fit prediction for the pure OPE calculation neglecting effects of the minimal photon energy 
cut on the OPE part (biases). The bands in figures (c) and (d) include those bias corrections of Ref. Q. The yellow bands 
represent the total experimental and theoretical fit uncertainty while the green band indicates the experimental uncertainty 
only. Solid markers are included in the fit while open markers are only overlaid for comparison. Moment measurements at 
different E cu t are highly correlated. 
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TABLE II: Results for the combined fit to all moments with experimental and theoretical uncertainties. For |V c t,| we add 
an additional theoretical error stemming from the uncertainty in the expansion for Fsl of 1.4%. Below the fit results the 
correlation matrix is shown. 



Combined 






OPE FIT RESULT: 


X 2 /N dof 


= 19.3/44 






Fit 


\V cb \ xlO" 3 


m b ( GcV) 


m c ( GcV) 


A (GeV 2 ) p% (GeV 3 ) p? G (GcV 2 ) p\ s (GcV 3 ) BR ct , (%) 


RESULT 


41.96 


4.590 


1.142 


0.401 


0.174 


0.297 


-0.183 


10.71 


A exp 


0.23 


0.025 


0.037 


0.019 


0.009 


0.024 


0.054 


0.10 


A HQE 


0.35 


0.030 


0.045 


0.035 


0.022 


0.046 


0.071 


0.08 


A F S L 


0.59 
















\Vcb\ 


1.000 


-0.399 


-0.220 


0.405 


0.267 


-0.305 


0.056 


0.700 


m b 




1.000 


0.951 


-0.387 


-0.189 


0.074 


-0.223 


0.098 








1.000 


-0.408 


-0.246 


-0.329 


-0.124 


0.143 


A 








1.000 


0.685 


0.257 


-0.008 


0.122 












1.000 


-0.050 


-0.479 


-0.055 


Mg 












1.000 


-0.035 


0.046 


Pis 














1.000 


-0.052 
1.000 




FIG. 4: Comparison of the different fit scenarios. Figure (a) shows the A\ 2 = 1 contour in the (rrib,fii) plane for the combined 
fit to all moments (solid red), the fit to hadron and lepton moments only (dashed blue) and the fit to photon moments only 
(dotted green). Figure (b) shows the results for the combined fit (solid red) and the fit to hadron and lepton moments only 
(dashed blue) in the (mt,,| V c {,|) plane. 
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IV. TRANSLATION OF FIT RESULTS INTO 
OTHER SCHEMES 

We translate the results for m b and /x 2 in the kinetic 
scheme to heavy quark distribution function parameters 
in other schemes so that they can be used for the ex- 
traction of \V ub \. The translation is done by predicting 
the first and second moment of the photon energy spec- 
trum above Scut = 1.6 GeV based on the heavy quark 
parameters from Table [H] and using the calculations of 
Rcf. 0. 

The experimental and theoretical uncertainties in the 
fitted parameters as well as their correlations are prop- 
agated into the errors on the moments as described in 
Section Till CI The minimum photon energy of 1.6 GeV 
is chosen such as to be insensitive to the distribution 
function itself. At this threshold the local OPE calcu- 
lation is applicable as the hardness Q = tub — 2£ , cut of 
the process is sufficiently high such that cut-induced per- 
turbative and non-perturbative corrections or biases are 
negligible. The predicted moments are given in Table ITTT1 

As the moments are physical observables which are 
scheme independent they can be used to extract the cor- 
responding heavy quark distribution function parameters 
in other schemes. For this translation, grids for the first 
and second moments of the photon energy spectrum are 
generated as a function of the two parameters (A,Ai) 
for Kagan-Neubert [3(| and (m b sFiP- 2 sf) f° r tnc Shape 
Function |33| scheme. A \ 2 is calculated for every set of 
parameters p = ((E^)(m b , /4), ((-E-y - (£ 7 }) 2 }(mf,, pi)) 
as 

X 2 = (Vi-Vi) V ij 1 (Vi-V-j) With V U = (J ^ a ]P^3 ( 7 ) 

where the yi are the predicted moments with their errors 
<7j and pij is the correlation between them. 

From the minimum value Xmin we obtain the central 
values for the parameters in the other schemes and de- 
termine the A\ 2 = 1 contour with respect to Xmin- 



TABLE III: First and second moment of the photon spectrum 
predicted for E cut = 1.6 GeV on the basis of the fit results for 
the HQE parameters. 

£ cut (GeV) <£ 7 }(GeV) {(E 1 — (£' 7 }) 2 }( GeV 2 ) p 
1.6 2.284 ± 0.018 0.0428 ± 0.0032 -0.03 



TABLE IV: Comparison of heavy quark distribution function 
parameters in the kinetic, Kagan-Neubert and Shape Func- 
tion scheme together with their correlation p. 

Kinetic Scheme 
7n b (GeV) M 2 (GeV 2 ) p 
4.590 ±0.039 0.401 ±0.040 -0.39 

Kagan-Neubert Scheme 
A(GcV) Ax (GeV 2 ) p 
0.621 ± 0.041 -0.497 ± -0.17 

Shape Function Scheme 
m bSF (GeV) /i 2 SF (GcV 2 ) p 
4.604 ± 0.038 0.189 ±0.038 -0.23 



Mathematica notebook based on Ref. EE Ell that 
was provided to us by the authors. In this calculation 
the moments are determined from a spectrum that is 
obtained by convoluting a shape function with a pcrtur- 
bative kernel with next-to-leading order accuracy, where 
we use the exponential form for the shape function given 
in Ref. |4jj. This calculation is conceptually similar to 
the one for B — * X u iv decays also presented in Ref. 0] 
which at present is used for the extraction of \V ub \ by 
several experiments. It therefore allows for a consistent 
determination of the shape function parameters for both, 
B — > X s j and B — ► X u iv decays. The numerical results 
for the shape function parameters are shown in Table llVl 
and the A% 2 = 1 contours are displayed in Figure 



A. Kagan-Neubert scheme 

In order to derive shape function parameters from the 
predicted moments in the Kagan-Neubert scheme |36j a 
grid of moments was generated for varying values of A 
and Ai. This was obtained using the Kagan-Neubert 
B — > A s 7 generator with the exponential Shape Func- 
tion ansatz as implemented in the Babar Monte Carlo 
generator (EvtGcn 38]). The results of this translation 
are shown in Table IIVI and Figure [SJ 



B. Shape-Function scheme 

For the translation into the Shape-Function 
scheme [3^ . l39j we use a grid of moments obtained with a 



V. APPLICATIONS FOR IMPROVED HEAVY 
QUARK PARAMETERS 

A. Improved OPE Expression for \V u b\ 

The results in the kinetic scheme for m b , // 2 , p? G and 
p 3 D have been used to give an updated expression for the 



standard local OPE formula for \V ub \ of Ref. mj,|4; 



\V ub 



4.268 • 10" 



j BR(B -> X u iv) l.&lps 



/ 0.002 
x(1±0.012 qcd ±0.022 H qe) 



tb 



(8) 



The error labeled 'QCD' includes perturbative uncertain- 
ties and those from weak annihilation. However, in con- 
trast to Ref. 0] we also include cxplicitely the 1 / con- 
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FIG. 5: Translation of fit results in the kinetic scheme to Kagan-Neubert (a) and the Shape Function scheme (b) via predicted 
photon moments. Figure (a) shows the results for the shape function parameters in the (A, Ai) plane from the combined fit to 
all moments (solid red) and the fit to hadron and lepton moments only (dashed blue). Figure (b) shows the corresponding fit 
results in the (mb sf,^ sf) plane. 



tribution from the p^ D term which results in an improved 
'QCD' error. The 'HQE' related uncertainty stems from 
the errors on mb, /x 2 ., an( ^ Pd anc ^ takes correlations 
between the parameters into account (see Tabic ITTf) . 

B. Extrapolation Factors for Measured B — > X s y 
Branching Fraction 



averaged where the total error was determined by com- 
bining the largest error from the scan of the error ellipses 
with half the maximum difference between any two mod- 
els in quadrature. Figure [5] also shows the spectra corre- 
sponding to the central values of Table llVl or cquivalcntly 
to the predicted photon energy moments of Table lllll in 
the three schemes. 



The measurement of the B — > X s j branching frac- 
tion is experimentally very challenging and has only been 
achieved for photon energies above -E cu t = 1.8 — 2.0 GcV. 
On the contrary, theoretical calculations predict the B — > 
X s j branching fraction at much lower values of E cut in 
order to avoid any dependence on the heavy quark distri- 
bution function. It is therefore customary to extrapolate 
measured branching fractions down to a value of 1.6 GeV 
where the y c an be compared to the theoretical calcula- 
tions 0, l46l | . Based on the heavy quark distribution 
function parameters in Table IIVI and the corresponding 
spectra we calculated a consistent set of extrapolation 
factors 



R(E cut ) — 



BR(B -> X sl ) Emt 

BR(B -> X,7)i.6GeV 



(9) 



for the kinetic, Kagan-Neubert and Shape-Function 
scheme. The results are summarised in Table IVl and Fig- 
ure The error was determined as the largest deviation 
from the central value obtained from a scan around the 
ellipses in Figures 01 and [5] where positive and negative 
errors were of comparable size. The results have been 



VI. CONCLUSION 

We have performed a fit to moments measurements 
from B — > XJLv and B — > X^y decays using calculations 
in the kinetic scheme S El The fit uses all currently 
available moment measurements from the BABAR Belle, 
CDF, CLEO and DELPHI experiments that are publicly 
available with their corresponding correlation matrices. 
We find that all the moment measurements of different 
order and from different inclusive B decays can be de- 
scribed by the fit result which is an important test of 
the consistency of this theoretical framework. We have 
extracted values for the CKM matrix element \V c b\, the 
quark masses nib and m c , and the kinetic expectation 
value /i 2 of 



\V cb \ = (41.96 ± 0.23 CX p ± 0.35 H qe ± 0.59 rsL ) x 10- 

m b = 4.590 ± 0.025 cxp ± 0.030 H qe GcV 

m c = 1.142 ±0.037 cxp ± 0.045 H qe GcV 

Hi = 0.401 ± 0.019 oxp ± 0.035 H qe GcV 2 
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FIG. 6: Figure (a) shows the photon energy spectra corresponding to the fitted heavy quark distribution parameters in the 
kinetic scheme (dashed red), Shape Function scheme (solid blue) and Kagan-Neubert scheme (dotted green). Figure (b) shows 
the corresponding extrapolation factors R(E cu t) for varying E cut for the kinetic scheme (red stars), Shape Function scheme 
(blue triangles) and Kagan-Neubert scheme (green squares) together with our average (black circles). 



TABLE V: Extrapolation factors R{E cu t) for BR(B -»• X s j). 





Kinetic 


Kagan-Neubert Shape Function 


Average 




Scheme 


Scheme 


Scheme 




£ cu t(GeV) 


R(Ecut) 


R(Ecut) 


^(-Ecut) 


R{E cu t) 


1.7 


0.986 ± 0.001 


0.988 ± 0.002 


0.982 ± 0.002 


0.985 ± 0.004 


1.8 


0.968 ± 0.002 


0.970 ± 0.005 


0.962 ± 0.004 


0.967 ± 0.006 


1.9 


0.939 ± 0.005 


0.940 ± 0.009 


0.930 ± 0.008 


0.936 ± 0.010 


2.0 


0.903 ± 0.009 


0.892 ± 0.014 


0.888 ± 0.014 


0.894 ± 0.016 



where the first error includes statistical and systematic 
experimental uncertainties and the second the theoreti- 
cal uncertainties from the HQEs. 

As can be seen, the error on \V c b\ which is below 2% 
is dominated by theorectical uncertainties. Any further 
improvements will require additional work on the accu- 
racy for the expression of Tsl, in particular on perturba- 
tive corrections to the Wilson coefficients of the the chro- 
momagnetic and Darwin operators. Similar observations 
can be made for mj,, which is determined with better than 
1% accuracy. However, the extraction of these quantities 
at the percent level represents in itself a remarkable test 
and success of the QCD-based calculations. 

The values for rrif, and /i 2 have been translated into the 
Kagan-Neubert Scheme where we obtain following values 
for the Shape Function parameters: 



A 

Ai 



0.621 ±0.041 GeV 
-0.497l°;°oSGeV 2 



Similarly, we obtain 



m b sf 

2 

Mtt SF 



4.604 ± 0.038 GeV 
0.189 ±0.038 GeV 2 



in the Shape Function scheme. As these parameters 
are critical for the extraction of \V u b\, their reduced 
uncertainty will enable measurements of \V u b\ at the 
5% level. This, together with \V c b\, will provide for a 
competitive measurement of the side of the Unitarity 
Triangle opposite the angle /3, and give further insights 
into the extent of CP violation in tree processes. 
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